Focusing of geodesic congruences in an accelerated expanding Universe 
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We study the accelerated expansion of the Universe through its consequences on a congruence of 
geodesies. We make use of the Raychaudhuri equation which describes the evolution of the expansion 
rate for a congruence of timelike or null geodesies. In particular, we focus on the space-time geometry 
contribution to this equation. By straightforward calculation from the metric of a Robertson- Walker 
cosmological model, it follows that in an accelerated expanding Universe the space-time contribution 
to the Raychaudhuri equation is positive for the fundamental congruence, favoring a non-focusing 
of the congruence of geodesies. However, the accelerated expansion of the present Universe does not 
imply a tendency of the fundamental congruence to diverge. It is shown that this is in fact the case 
for certain congruences of timelike geodesies without vorticity. Therefore, the focusing of geodesies 
remains feasible in an accelerated expanding Universe. Furthermore, a negative contribution to the 
Raychaudhuri equation from space-time geometry which is usually interpreted as the manifestation 
of the attractive character of gravity is restored in an accelerated expanding Robertson- Walker 
space-time at high speeds. 

PACS numbers: 04.50.-h, 98.80.-k, 04.20.-q 



I. INTRODUCTION 

General Relativity (GR) as a theory that describes the 
geometrical structure of space-time and its connection 
with the energy content of the Universe has been highly 
successful in the last one hundred years. GR is not only 
able to give explanation of some discrepancies that had 
been arisen with Newtonian theory in the Solar System, 
but also provides a satisfactory description of the cosmo- 
logical evolution of the space-time. 

However, GR is not able to account for the latest 
cosmological and astrophysical observations with stan- 
dard matter sources. Among these observations are the 
anomalous rotational velocities of objects near the edge 
of the galaxies, the dynamics of galaxies in clusters, for- 
mation of large scale structures or gravitational lensing 
of background objects by galaxy clusters, such as the 
Bullet Cluster. For the explanation of these observa- 
tions, more matter than the standard one is required, 
the so-called dark matter (DM) which has an attractive 
gravitational character that decelerates the expansion of 
the Universe. Although there are many plausible origins 
for this component DM is usually assumed to be in 
the form of thermal relics that naturally freeze-out with 
the right abundance in many extensions of the standard 
model of particles @] • Future experiments will be able to 
discriminate among the large number of candidates and 
models, such as direct and indirect detection designed 
explicitly for their search or even at high energy col- 
liders, where they could be produced Furthermore, 
the authors of [B[ , recently awarded with the Nobel Prize 
in physics, discovered the accelerated expansion of the 



Universe. This result has been one of the main reasons 
for the revision of the cosmological evolution as predicted 
by GR. Neither the properties of usual matter and en- 
ergy nor the believed properties of DM can account for 
this acceleration. One way of circumventing this prob- 
lem is to postulating a new kind of fluid, the so-called 
dark energy, which has a non-standard equation of state 
allowing a positive energy density while having a nega- 
tive pressure 0, 0| ■ Another possibility to generate the 
present accelerated expansion of the Universe is by the 
modification of GR [8|, [9( . Indeed, this type of modifi- 
cations could produce an inflationary epoch in the very 
early Universe. This epoch is postulated to evade the 
problems with the initially homogeneous and isotropic 
state of the Universe which arise from considering parti- 
cle horizons in Robertson- Walker cosmological models. 

In this work, we will analyze the accelerated expan- 
sion of the Universe by its description through its effects 
on a congruence of geodesies without assuming any spe- 
cific gravitational theory. This description has the ad- 
vantage of being coordinate-independent. In an acceler- 
ated expansion it is expected that neighboring geodesies 
with similar direction will increasingly recede from each 
other. This evolution requires the definition of a kine- 
matic quantity, the expansion, which is quite analogous 
to the spatial divergence of a vector field. The expan- 
sion rate along the congruence will provide us with the 
appropriate tool to analyze the consequences of the cos- 
mological accelerated expansion of the Universe. This 
rate is given by the Raychaudhuri equation (RE), which 
will be discussed and used throughout this work. RE was 
first derived by Raychaudhuri in [lpl ] for congruences of 



timelike geodesies. Later, Sachs introduced the cor- 
responding RE for null geodesies. Ehlers et al. [IH pro- 
vided a generalization of hydrodynamics to GR where re- 
sults analogous to well-known theorems of hydrodynam- 
ics using RE equation were obtained. Given a congruence 
of timelike or null geodesies this equation provides the 
contraction or expansion rate of the congruence accord- 
ing to some kinematic properties of the congruence and 
the curvature of the space-time where the congruence is 
embedded. In GR, the RE together with assumed energy 
conditions guarantee the attractive character of gravity 
represented by a non-positive contribution to this equa- 
tion from the space-time geometry. Furthermore, the RE 
has a fundamental role in the demonstration of important 
theorems due to Hawking and Penrose. These theorems 
assert that singularities of manifolds representing phys- 
ical space-times are intrinsic to GR and not a special 
property of particular solutions of Einstein's equations 
provided that some energy conditions hold HI, T3 |. 



In this investigation, we shall analyze the evolution of 
the expansion of a congruence of geodesies. This is car- 
ried out in a Robertson- Walker cosmological model from 
a kinematical point of view since no gravitational theory 
is assumed beyond the equivalence principle. This work 
is organized as follows: First, in Section [TTl we present a 
short review of the derivation of the RE for illustrative 
purposes and stress some aspects usually disregarded in 
the available literature. Then, the kinematics of a con- 
gruence of geodesies in a Robertson- Walker space-time is 
studied in Section iHll We conclude the paper by giving 
our conclusions in Section IV. 



Congruences of timelike geodesies 

Let us consider a congruence of timelike geodesies and 
its tangent vector field £ a . It is assumed that the tan- 
gent vector field is normalized, i.e., £ a £ a = — 1, which is 
always feasible by choosing an affine parameterization of 
the congruence. Let us define the (0,2) tensor 



^ab = Vb^ a 



It is easy to verify that 
C^ab = o 



^ab = 0. 



(1) 



(2) 



The first expression is a consequence of £ a being normal- 
ized whereas the second one is the geodesic equation. As 
usual the projector tensor h a b is defined by 



hab — 9ab + £,a£,b ■ 



(3) 



This is a spatial metric provided that £ a is hypersurface 
orthogonal. This property of the vector field means that 
there exists a foliation of the manifold by hypersurfaccs 
orthogonal to £ a . The metric induced by g a b on these 
hypersurfaces is h ao . For cases where an orthogonal hy- 
persurface cannot be found, the projector tensor h a b is 
still well-defined as the tensor that projects tensors in 
a subspace orthogonal to £ a at each point of the mani- 
fold. It is worth noticing that S a b has no components in 
the direction of £ a according to expression (|2|). Conse- 
quently, S a b is rather defined in the subspace orthogonal 
to £ a . This subspace is the hypersurfaces tangent space 
foliating the manifold when £ a is hypersurface orthogo- 
nal. Finally, let us consider the following quantities 



~ ab h 

i — i I In 



(4) 



II. RAYCHAUDHURI EQUATION 
DERIVATION 



&ab = S( ab ) - - 9h ab , 



(5) 



The RE provides the expansion rate of a congruence 
of timelike or null geodesies. In the following, we will 
summarize the main results for both cases separately ac- 
cording to the notation in 14( . For different reviews and 
approaches to the RE, we refer the reader to [Hj]. In 
the following let us use abstract index notation (see [14 1 
for details). Thus, an expression like £ a represents the 
vector itself and not a component in a certain basis. An 
expression like g ab £, a (j) b must be understood as g (£,</>)■ 
When explicit mention to the components of a tensor in 
a certain basis is made, greek indices will be used, i.e., 
The principal advantage of this notation is that it 
distinguishes equations that hold between tensors and 
equations that only hold for components in a certain ba- 
sis. We use a metric signature (—,+,+,+) and the Rie- 
mann and Ricci tensors and the Ricci scalar are defined 
respectively as R abc d = d b T d ac - d a T d ( -J- re rd - re rd 
R ac = R a0 c and R = R a a - 



Uab = ^[ab] 



(6) 
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which are the trace (with respect to the metric h a b, but 
the trace is the same respect to g a b because of ([!])), the 
trace-free symmetric part and the antisymmetric part of 
S Q b respectively. These quantities are known as the ex- 
pansion 9, shear a a b and twist u ab of the congruence. 
Thus, S a b can be decomposed as 



2<j& — - Oh a b + dab 



■ W a b , 



(7) 



which is analogous to the decomposition of the deforma- 
tion velocity tensor in hydrodynamics. The meaning of 
9, cr a b and u; a b is clear from the analogy with hydrody- 
namics. The expansion 9 takes into account the diver- 
gence or convergence of nearby geodesies. It may also 
be interpreted as a relative change in volume of a fluid 
particle. The shear er a b provides information about the 
deformation of the geodesic congruence without change 



2 



in volume. Finally, the twist uj a b contains the vorticity 
of the congruence. Equations for the rate of these kine- 
matic quantities along each geodesic can be derived in 
the following way. After some calculus, it can be proved 



(8) 



where the Leibniz's rule and the geodesic equation satis- 
fied by £ a have been used. Since we are interested in the 
equation for the expansion 6, the trace of (j8]) becomes 



d6 
d7 



^ = - o ah G ab + uj ab u ab - R ab ee 



(9) 



where t is an affine parameter of the geodesies, i.e., the 
proper time of an observer moving along a geodesic of the 
congruence. Equation © is known as the RE [ToL [lH . 
Let us recall that the expansion provides a natural 
coordinate-independent way of measuring the focusing or 
divergence of free-falling matter following a certain con- 
gruence of geodesies in the gravitational field, i.e., 6 > 
means divergence and 6 < means convergence. The 
quantities of the r.h.s. of (|9]) are scalars, thus any conve- 
nient coordinate system can be used to calculate them. 
However, it is important to stress that equation © pro- 
vides the expansion rate as seen by an observer moving 
along a geodesic of the congruence. The equivalent coun- 
terparts for shear a ab and the twist io a b are obtained from 
equation by taking the trace-free symmetric part and 
the antisymmetric part respectively. 

Let us now discuss the character of the contribution 
of the different terms in expression ([9]) . The quantities 
Cabf afc and u! a i,uj ab are non-negative quadratic scalars. 
The fact that these scalars are non-negative is due to the 
spatial character of the tensor S a & reflected by ([2]). This 
can be understood if one considers a framework with the 
spatial components of h a b, i.e., the orthogonal compo- 
nents to £ a , taking the values 8^ v with /x, v = 1,2,3 at 
a certain point of the manifold, whereas the components 
in the direction of £ a being zero. Then, taking into ac- 
count the spatial character of S a b the values of a^v and 
are the same in this special framework. Therefore, 
PabO' ab > 0. The same reasoning is valid for uj a b- 

Thus, it is clear that both the shear a a b and the ex- 
pansion 6 always tend the congruence to focus whereas 
the twist LUab always tends to make it diverge. For the 
contribution of space-time geometry, i.e., for the sign of 
—Rab£, a £, b , there appear three possibilities: 

- Rabd b > Positive contribution, (10) 



Rabd b < Negative contribution, (11) 



- RabCt = Zer0 contribution. (12) 

The implications of the last three inequalities will be the 
keystone of this investigation. The effect that the two 



first possibilities produce in the congruence may be ex- 
pressed as a tendency to diverge and a tendency to focus 
respectively. It can also be understood as accelerated ex- 
pansion in the case of (|10j) and decelerated expansion in 
the case of (TTTI) . On the other hand, the value of i? a &£ a £ a 
depends both on the point of the manifold where it is 
calculated and on the direction of the vector field £ a at 
that point. Thus, a positive contribution to the RE may 
be obtained at some points of the manifold for certain 
directions whereas may be negative at those points for 
different directions. 

Although it is rarely considered, the expression 
Rab£, a t; b admits a geometrical interpretation [16|. Let 
£?q% be a unit vector and £?s with a = 1,2,3 an arbi- 
trary set of unit vectors orthogonal to £" j . The Gaussian 
curvature K(o<x) 01 t ne geodesic surface generated by £" ) 
and £ ( a Q) is 

K, 

where e(o) 

take the values +1 or — 1 according to ffafctfo) £(b) = e (°) 
and g ab £ ( q q) = e (ct) . As is widely known, the sign 
convention for Riemann and Ricci tensors is arbitrary a 
priori. However, let us recall that the sign of the Gaussian 
curvature possesses an intrinsic meaning sign-convention 
independent. Thus, when all these Gaussian curvatures 
K(n a -\ are summed one gets 



(On 



e(0) e( a ) Rabcd £(0) £(a) £(0) £(at) ' ( 13 ) 



and 



<<*) 



are the usual Eisenhart e's which 
-1 according to g a b^ 0) £ ( b 0) 



a=0 



K, 



(Oct) — e(Q) Rac£(p) £(0)' 



(14) 



where we have considered that if(oo) is zero and the spec- 
tral decomposition of the metric 



ab 



3 

E 

a=0 



e (a) i( a ) 



(15) 



has been used. Manifestly, the sum (114|) does not depend 
on the election of the set of orthogonal unit vectors £° a j 
under consideration. Then, eR a b£, a £, b is the sum of all the 
Gaussian curvatures of the geodesies surfaces generated 
by £ a and any set of orthogonal unit vectors to it. This 
has been referred to as the mean curvature in the direc- 
tion of f a (l6l]. Thus, by following our sign convention 
e = — 1 and the mean curvature in the direction of 
a timelike vector £ a becomes 

= -R ab e^ ■ (16) 

Hence, the contribution of the space-time geometry to the 
RE ([9]) has a clear geometrical interpretation as the mean 
curvature Ai^ in the direction of the tangent vector field 
£ a . Furthermore, the inequalities (fT0|) . (JTTJ) and (fT"2"|) can 
be rewritten in terms of the mean curvature M^a in the 
direction of £ a as follows 



Mt* > 



Positive contribution. 



(17) 
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M.^a < Negative contribution, 



Me 







Zero contribution. 



(18) 



(19) 



Thus, if the mean curvature for a given timelike vector 
£ a at one point of the manifold is positive (negative) the 
space-time geometry would tend to make geodesies with 
tangent vector £ a passing through this point diverge (fo- 
cus). This formulation represents a clear advantage over 
the previous one in terms of R a b£, a £, a since it does not 
depend upon any sign convention nor upon the metric 
signature. Let us recall that when energy conditions are 
assumed in GR, one always obtains (fT8| and (fT9]) but 
not (fT"7|) . The fact that in GR is always negative 

or zero is interpreted as a manifestation of the attractive 
character of gravity. 

Let us return to our discussion about the RE. It is 
clear that the sign of the expansion rate dO/Ar will de- 
pend upon all the contributions of the r.h.s of (|9|). From 
their respectively definitions, it is obvious that all these 
terms depend on the congruence of geodesies under con- 
sideration. However, the terms 8 2 , a a b(J ab and u> a b^ ah 
depend on the behavior of the congruence in the neigh- 
borhood of the point where they are calculated. Let us 
recall that these terms come from V{,£ a , then the tangent 
vector field £ a must be known locally in order to calculate 
its covariant derivative. Nevertheless, the contribution of 
space-time geometry depends only upon the direc- 
tion of the congruence at a given point. Thus, all the 
congruences pointing in the same direction at a given 
point, will have the same contribution from space-time. 
Moreover, there are only two terms that may have a pos- 
itive contribution: the term oj a bUJ ab and the one given 
by M^a. Specifically, we shall restrict ourselves to con- 
gruences with uj a b = as we shall explain at the end of 
this section. Therefore, in this type of congruences the 
only term that may have a positive contribution to (|9]) is 



Congruences of null geodesies 

An equivalent version of the RE (|9]) for null geodesies 
can be obtained by considering the null tangent vector 



field k a to a congruence of null geodesies 
expected the result is 



lMl. As 



d0 
dA 



~ Z-ab i ~ ,~,ab 



R ab k a k b , (20) 



where the quantities a a b, &>ab and 9 are defined with re- 
spect to the congruence of null geodesies, the tensor h a b 
is the tensor analogous to h a b for null geodesies and A 
is an affine parameter in the following sense. When null 
geodesies are considered the tangent vector field has al- 
ways zero norm and consequently an affine parameter 



defined as the one which normalizes the tangent vec- 
tor field does not hold. For this reason, an affine pa- 
rameter can be defined in a more general way as the 
one preserving the geodesic equation in its usual form. 
This definition remains valid both for timelike and null 
geodesies. The factor 1/2 in (|2"0|) instead of 1/3 appear- 
ing in ([9]) is due to the fact that the involved subspace in 
the case of null geodesies is two-dimensional rather than 
three-dimensional. Thus, equation (|20|) is the RE for null 
geodesies. 

Finally, let us remind that the same type of space- 
time geometry contributions provided by the inequali- 
ties (HHD, dTTJ and (H2) applies replacing R a b£, a ^ b by 
R a bk a k b . However, an analogous interpretation of the 
term R a bk a k b as a mean curvature is not practicable in 
this case 11611. 



Further remarks 

At this stage, we estimate important to stress that the 
Raychaudhuri equations ([9]) and (|20[) are identities, i.e., 
geometrical statements which do not require any under- 
lying gravitational theory to be derived. They are quite 
general and apply whenever a RiemanniarQ geometry for 
the space-time is assumed. 

A more general equation for the evolution of the ex- 
pansion of a congruence results if the assumption of being 
a congruence of geodesies is discarded. Moreover, the RE 
can be generalized to spaces with torsion 17[. For our 
purposes ((9|) and (|20|) are sufficient. 

As was briefly mentioned, there are analogous equa- 
tions to ([9]) and ([20| for the shear a a b{a a b) and twist 
oj a b(u>ab) that provide the rate of these quantities along 
geodesies. These equations form a coupled system of non- 
linear first order differential equations. They also receive 
the generic name of Raychaudhuri equations. For in- 
stance, the equation for the twist ui a b is an homogeneous 
differential equation in to a f, and thus if uj a b is zero initially 
it will remain zero along the geodesies. This result does 
not prove that the unique solution is u> a b = because of 
the non-linearity of the equation. However, this conclu- 
sion can be eventually proved more rigorously (see [12J ) . 
This is in fact the analogous result to Kelvin's theorem 
in hydrodynamics showing the analogy between theorems 
of hydrodynamics and theorems concerning congruences 
of geodesies [ijj . 



The initial value problem of Raychaudhuri equations is 
still an open question. However, in the case of equation 
© if congruences that satisfy uj a b = are considered, 



1 Let us stress that we use the adjective Riemannian in an extended 
sense, i.e., Lorcntzian metrics are included. 
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there would exist hypersurfaces orthogonal to the vector 
field that foliate the space-time. This is a consequence 
of Frobenius' theorem [13]. Therefore, for a normalized 
vector field defined orthogonally to an arbitrary hyper- 
surface - not necessarily in the whole manifold - its shear 
and its expansion (with components only in that hyper- 
surface ([2])) can be determined. This way, the expansion 
rate of the congruence of geodesies defined by this vector 
field can also be determined. 

Therefore, the rate of the kinematic quantities of the 
congruence (9, Lo a b and o~ab for timelike geodesies and 9, 
ujab and a a b for null geodesies) is influenced by the geome- 
try of space-time. In particular, in the case of expansion 
9 (6), the term -R a bC£ b (-R a bk a k b ) is the contribu- 
tion of space-time geometry for timelike geodesies (null 
geodesies). Let us recall that the contribution of space- 
time geometry has more general implications than the 
other contributions since the former depends solely on 
the direction of the tangent vector field at one point but 
not on the local behavior of the tangent vector field as 9 
(9), a a b (o'ab) and U! a b (&ab)- Thus, the contribution of 
space-time at one point of the manifold is the same for 
every congruence of geodesic with the same value of the 
tangent vector field at that point. However, the quanti- 
ties 9, a a b and uj a b will be different if the congruences are 
not locally the same. 

Finally, let us remind here that the field equations de- 
rived from the assumed gravitational theory determine 
the geometry of the space-time. For instance, in the case 
of GR, the usual Einstein's field equations enable to re- 
late the quantities R a b^ a ^ b and R a bk a k b which appear in 
([9]) and (|20j) respectively with the energy-momentum ten- 
sor T a b in a very simple way. Hence, the expansion rate 
is straightforwardly related with the energy-momentum 
content. As a result, general properties of the behavior of 
the geodesic congruence can be established provided that 
some energy conditions on the energy-momentum tensor 
are assumed Lj, LL4[ . More precisely, if the strong energy 
condition and the null energy condition are assumed (al- 
though the last one also follows by continuity if the strong 
energy condition holds) one reaches the geodesic focusing 
theorem, which assures the convergence until zero size of 
a congruence of geodesies without vorticity which is ini- 
tially focusing. This result is used in the demonstration 
of the singularity theorems by Hawking and Penrose and 
it has recently drawn wide attention due to its prominent 
role in the holographic principle (l8| . 



III. KINEMATIC CONSIDERATIONS IN 
ROBERTSON- WALKER SPACE-TIMES 



In this section, let us start by considering a spatially 
flat Robertson- Walker metric [19(, 



that according to the ultimate measurements of cosmo- 
logical densities of the Universe [2(J , represents correctly 
the Universe at sufficiently large scales in agreement with 
the hypotheses of homogeneity and isotropy stated in the 
Cosmological principle. The connection between the in- 
equality (TlU)) - or equivalently <jTTJ) - with the accelerated 
expansion of the Universe is not straightforward. The 
accelerated expansion of the Universe in a Robertson- 
Walker cosmological model is governed by the condition 
a > where the derivatives are taken with respect to 
cosmic time. However, provided that the inequality (|17[) 
holds for all the timelike directions £ a in a region of 
the manifold, the contribution of space-time geometry 
then favors a tendency to diverge of all the congruence of 
geodesies in that region. Therefore, the natural question 
that arises is whether the condition a > 0, when satisfied 
in a region, implies M.^a > for every timelike direc- 
tion £ a . This question will be addressed throughout the 
present section. 

One may thus wonder about the effects of the cos- 
mological acceleration in the kinematical behavior of a 
congruence of timelike geodesies. In particular, it is of 
special interest to study whether a congruence of time- 
like geodesies tends to focus - or to diverge - and under 
which conditions this happens. The kinematic quantity 
of a congruence known as expansion 9 introduced in Sec- 
tionjniprovides us with a useful tool to study the focusing 
of timclikc geodesies. In this sense, the inequality 

c\9 

£°V C = -=- > , (22) 

QT 

states the tendency of nearby geodesies to diverge mea- 
sured by a comoving observer. The proper time r is the 
time of the comoving observer. Since the l.h.s. of the pre- 
vious expression (|22l) can be replaced by the r.h.s. of the 
RE ©, a vector field tangent to a congruence of time- 
like geodesies must be chosen in order to determine the 
expansion 9, shear a a b and twist ui a b- 



Fundamental congruence 

Since we are dealing with a Robertson- Walker metric 
211) , the most natural choice for a vector field would be 



few - dt , 



(23) 



which is the preferred direction of a Robertson- Walker 
space-time because of the local isotropy in hypersur- 
faces orthogonal to fe w . The vector field (|23| gives 
rise to a congruence of timelike geodesies of the form 
x^ = (t — to, xq) in comoving coordinates with to and xq 



constants [ff These geodesies are followed by the funda- 



ds 2 



-dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) 



(21) 



2 Therefore, being steady is a geodesic motion in Robertson- Walker 
metric. 
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mental observers, those seeing homogeneity and isotropy 
in constant-time hypersurfaces. Thus, a fundamental 
congruence refers to the congruence of timelike geodesies 
generated by ^ w . As it is easy to verify by straightfor- 
ward calculation - or by using symmetry arguments -, 
the kinematic quantities of the fundamental congruence 
are all zero except the expansion 



'RW 



3^. 

a 



(24) 



Thus, by inspection of © it follows that the rate of the 
expansion seen by a fundamental observer will be entirely 
determined by the value of the expansion #kvv and the 
contribution of space-time geometry 



Hence 



d6 



RW 



df 



RW 



a \ a 
- +- 
a I a 



(25) 



(26) 



where the proper time of the fundamental observer is t, 
i.e., the time coordinate in (|2~T1) . A fundamental observer 
will then see nearby geodesies with a tendency to diverge 
if and only if d#Rw/dt > 0. By imposing this condition 
on (j2"6"j) we get 



a > 



which can be expressed in the following way 



1 a 



< -1 



dd 



mv 



dt 



> 0. 



(27) 



(28) 



where q is the deceleration parameter and H = a/ a is the 
Hubble parameter. Hence, a fundamental observer will 
perceive nearby geodesies with a tendency to diverge if 
and only if q < — 1. In other words, getting a decel- 
eration parameter smaller than zero is not sufficient for 
fundamental timelike geodesies to be perceived by a fun- 
damental observer as diverging in an increasingly fashion. 

In the standard cosmological model ACDM, the value 
of q today can be determined by using the values of the 
present relative abundances of matter (including dark 
matter) and cosmological constant, £Im and Q\ respec- 
tively by the expression 



(29) 



where flatness in spatial sections was again assumed. 
Thus, using the approximate values flj^ ~ 0.27 and 
SI A ~ 0.73 given in [20(, the present value of the decelera- 
tion parameter go results go — —0.6, with an uncertainty 
of the order of 0.01. This value of go does not imply a ten- 
dency of nearby geodesies to diverge since g > —1- One 
may argue that we have assumed GR to calculate the 
deceleration parameter qo, and therefore the discussion 



is not independent of a gravitational theory. However, 
there are straightforward estimations of the deceleration 
parameter using Supernovae la [21| without appealing to 
a gravitational theory. The preferred values obtained for 
<7o in these analyses are always larger than —1, although 
with uncertainties of the order 0.1. 

Before continuing with the study of the kinematics of 
a congruence of timelike geodesies, let us mention that 
in general, the determination of other geodesies is not 
possible in a closed form since an expression for a(t) is 
usually missing. Therefore, it will not be possible in gen- 
eral to obtain an expression for the tangent vector field 
£ a of a congruence of geodesies. However, if our interest 
is focused on calculating the space-time geometry contri- 
bution to the RE for timelike geodesies, the knowledge of 
the vector field £ a expression is not necessary. The reason 
is that the value of only depends on the direction 
of the vector field at a given point. Then, given a time- 
like vector, the value of Ai can be determined by using 
the expressions for the components of the Ricci tensor. 
Moreover, the value of Ai^ will be constant in each hy- 
persurface orthogonal to £g_ w because of homogeneity in 
these hypersurfaces. 

However, if we wish to know the total value of d8/dr 
from the RE, an expression for the tangent vector field 
is needed. It means that we must know the vector field 
locally at each point of the manifold. Nevertheless, we 
can use the property stated in equations ([5]) which as- 
serts the spatial character of S a b = Vb£ a . The procedure 
consists on defining a smooth vector field on a spacelike 
hypersurface. The only conditions on the definition of 
this vector field are that it be normalized £ a £ Q = — 1 
and that it be orthogonal to the hypersurface where it is 
defined. This slice of a vector field will give rise to a con- 
gruence of timelike geodesies passing by the hypersurface 
with the direction of the vector field in each point. This 
congruence will have uj a b — because being zero in the 
initial hypersurface it will remain so. However, the shear 
Gab and the expansion 9 do not need to be zero. Then, 
instead of considering a vector field tangent to a con- 
gruence of geodesies defined in all the manifold, a vector 
field is instead defined at a spacelike hypersurface. Let us 
remark that this procedure restricts us to consider con- 
gruences of timelike geodesies without vorticity. By this 
way S a (, at the initial hypersurface can be determined. 
Thus, the expansion rate dO/dr of the congruence eval- 
uated in the hypersurface will be known from a slice of 
the geodesic tangent vector field. 



One-parameter dependent congruences of timelike 
geodesies 

The congruences under consideration from now on are 
obtained from a one-parameter dependent slice of a vec- 
tor field. First of all, we must choose a spacelike hy- 
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persurface where the vector field will be defined. The 
simplest one-parameter dependent e family of spacelike 
hypersurfaces defined by a function / = in terms of the 
coordinates (I2T1) is 



f = t-e-x = 0, 



(30) 



where < \e\/a < 1 in order the hypersurfaces to be 
spacelike. This constraint is easily understood if one con- 
siders the normal vector to the hypersurfaces (j3"0)) and 
this vector is then constrained to be timelike. This fam- 
ily of hypersurfaces depends in fact on three parameters, 
namely the three components of e. However, thanks to 
the spatial isotropy, the direction of e does not affect the 
results. Thus, the relevant parameter is |e*| = e. More- 
over, because of spatial homogeneity and the arbitrari- 
ness in the election of the time origin we have eliminated 
the constant that may appear in ([30|) . i.e., the hypersur- 
faces pass through the origin of coordinates. The above 
restriction on the value of e means that the hypersurfacc 
is defined only in the manifold region a > e. Otherwise, 
the hypersurfaces are not spacelike. Let us now define 
the following quantity 



a?{t) 



(31) 



The slice of the timelike vector field which is normalized 
and orthogonal to the hypersurface given by (|30l) then 
becomes 



ft 



(0) 
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d_ 

dt 



P 



d_ 

dx 



where 
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y/l -a?P 2 



(32) 



(33) 



The tangent vector field that gives rise to the fundamen- 
tal congruence is obtained in the specific case e = or 
equivalently j3 = 0. In order to guarantee to be time- 
like, one imposes < af3 < 1. Then, it follows the re- 
striction on the value of e commented above. From (1321) . 
/3 can be interpreted as the spatial coordinate velocity of 
a particle passing through the hypersurfaces defined by 
(f30|) . The dependence of (3 on the scale factor (|3Tj) means 
that this velocity is not the same in the whole hypersur- 
face defined by a certain value of the parameter e. Note 
that in expression Q32p we have chosen the cc-axis in the 
direction of the spatial coordinate velocity (3 without loss 
of generality thanks to the spatial isotropy. The 7 factor 
just ensures the normalization <? a 6£?g)£('s) = — 1- 

One realizes of the similarity of the former derivation 
with the usual definitions of f3 and 7 in Special Rela- 
tivity. Consequently, the slice of the timelike vector field 
orthogonal to the hypersurfaces (|30|) in Special Relativity 
becomes 



(34) 



with < (i < 1 in order to guarantee rj a bU a u b < 0, with 
r\ a b the Minkowskian metric. Hence, expressions (1321) and 
(f3"3"|) are completely valid both in Special Relativity (set- 
ting a — 1) and in Robertson- Walker space-times. How- 
ever, when we are considering a Robertson- Walker met- 
ric, the maximum coordinate velocity (3 max that a phys- 
ical entity can reach, in the sense that it moves along 
null geodesies, depends on the value of scale factor a(t) 
at each point. In order to guarantee that the vector field 
as given by (|32p is timelike, the required condition reads 
< aft < 1 . Therefore the maximum coordinate velocity 
is expressed as f3 max = l/a(t). Nevertheless, let us recall 
that the physical spatial velocity which is measured is 
a/3. One may interpret a physical velocity as the ratio 
between a spatial and time distances. These distances 
hold for the ones induced by the timelike geodesic one 
follows. Therefore, for a particle with four- velocity 
its physical velocity yields 



V, 



physical 



lim 

At->0 



rt + At 



CRWta Jj 
?a ?(/3) al 



a/3 (35) 



where w is given by flM|) and = g ab + £* W 4 RW 
is the spatial metric induced by the flow of time in the 
direction of . Alternatively, one may also build a unit 
vector orthogonal to being the spatial projection of 

^ namely h a = h™^/^^™^™^, and 
then calculate the ratio between the projection of in 
the direction of h a and in the direction of £r W , finally 
obtaining the same result as in (1351) . 

Therefore, after considering the slice of the timclikc 
vector field (|52"j) defined in the hypersurface (|5U)) . 
the tensor S a f, = Vt,£ Q is calculated in order to obtain 
both the expansion 9 and the shear a a b (the twist uj a b is 
zero by construction). By following this procedure, the 
value of d#/dr in the hypersurfaces / = t — e '■ x will be 
determined. This will provide us with information about 
the focusing or divergence of the congruence generated 
by the vector field slice in terms of the scale factor a(t), 
its derivatives and the spatial coordinate velocity /3. 

However, in the following we shall focus on the mean 
curvature Ai^a. This is due to the fact that although 
(Jab and 6 may not be zero, their contribution to the RE 
is always non-positive. Moreover, as stressed in Section 
im the implications of the mean curvature Ai^a are more 
general than those of the other terms. Furthermore, M.^ 
is the unique term that may have a positive contribution 
to the RE ([9]) since for the congruence given rise by the 
vector field slice ([32]) we have u) a buj ab — 0. Hence, let us 
proceed with the calculation of in order to extract 

information about the focusing of these one-parameter 
dependent congruences. 
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Evaluation of the condition Me a > 

In the following we will obtain a condition on a(t) 
and the spatial coordinate velocity /? that guarantees 

> ^ ' *' e '> a P os itive contribution to the RE (fT7| . 
The motivation of choosing a positive contribution lies 
in the fact that the fundamental congruence satisfies 

> in the present epoch as follows from (|25|) . 

Then, by using the Ricci tensor components derived 
from the metric (j2"Tj) and the expression for the timelike 
vector field slice given in (j52")l . we get 



-3- + [3 2 (ad + 2d 2 

a v 



(36) 



Thus, imposing Mq m = -Rab^S,^ > 0, yields 



3-+(3 2 (ad + 2d 2 ) <0 
a v 7 



j^d 2 



where a 2 /? 2 — 3 < (because of a/3 < 1) have been used. 
This last inequality tells us that the condition a > is 
not sufficient for a timelike vector of the form (|52")l to get 
a positive contribution to the RE from the space-time 
geometry M.^ . Consequently, d must be larger than a 
certain positive value. Therefore, a positive contribution 
to the RE from space-time geometry (jTTJ) for all timelike 
congruences is not equivalent to d > 0. The inequality 
(l37t can be recast as 



Q < 



-2a 2 /? 2 
3 - a 2 /3 2 



(38) 



which depends only on the physical velocity a/3. In Figure 
[U the condition (|38| has been depicted. In that figure, 
the values (a/3, q) ensuring a positive contribution to the 
RE from the mean curvature are plotted. 

Moreover, it is worthwhile to stress that the require- 
ment M.ea > implies lower values of the deceleration 
parameter q for timelike vectors approaching the null 
cone. Then, in the limit a/3 — > 1, i.e., timelike vectors 
near the null cone, we get q < — 1. Therefore, provided 
that q < —1 all the congruences of timelike geodesies will 
have a positive contribution to the RE from space-time 
geometry which means -M^ > for all timelike vectors 
£?gv The implication of this condition must not be con- 
fused with (J55]) : the condition (J35J) implied d6> RW /dt > 
for the fundamental congruence whereas in this case the 
condition q < — 1 implies a positive contribution from 
the space-time geometry, namely M-^a > 0, for every 
congruence of timelike geodesies. Since the present value 
of the deceleration parameter is not smaller than —1, it 
remains possible to get M.^ < for timelike congru- 
ences whose tangent vector field is sufficiently inclined 
with respect to £j| w = d/dt. Therefore, using (1551) . the 
following possibilities for the physical velocity a/3 may be 
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Figure 1. Mean curvature sign for congruence of geodesies gener- 
ated by i ] J >2[ ): The line represents the limiting case q = — 2a 2 /3 2 /(3 — 
a 2 /3 2 ). For a given spatial coordinate velocity /3 and scale factor 
a(t) if the deceleration parameter q is below the line (region in 
blue) a congruence of geodesies with a physical velocity a/3 ac- 
quires a positive contribution to the RE from the mean curvature 
M^a or equivalently from the term —R a b£, a £, b ■ lu the region above 
the line (white color) we depict the opposite case, i.e., negative 
contribution. 



fulfilled 



a/3 < 



a/3 > 



a/3 = 



-3g 
2-q 



Positive contribution, 



(39) 



-3q 



-3q_ 

- q 



Negative contribution, (40) 



Zero contribution; 



(41) 



to the RE from the space-time geometry. 

If j3 = in (f3"T|) we obtain a > which means that 
the fundamental congruence will experience a positive 
contribution from the space-time if d > and therefore 
in agreement with (|25p . 

In Figure [5J a future light cone is represented. Using 
the value for the present deceleration parameter (qo — 
—0.6), the limiting value obtained from (|4*T|) . a/3 ~ 0.8 
separates the regions of different contribution to the RE. 
The requirement a/3 > 0.8 constitutes a necessary con- 
dition to ensure M^a < 0. However, due to the non- 
positive character of all the other terms of the RE for 
a congruence without vorticity, dO/dr might be negative 
despite a positive contribution A4^a > 0. 
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Figure 2. Modifed future light cone in Robertson- Walker metric: 
The line t = ax corresponds to the usual light-cone. A positive 
contribution to the RE (i.e., M^a > 0), for timelike vectors £ a lies 
in the inner (red) area and expression H39II is there satisfied. For 
vectors in the outer (blue) zone M^a < and condition 1 140 II is now 
obeyed. The value of the slope of the straight line that separates the 
two regions corresponds to 1/(0,0) = 1/0.8. This value is obtained 
using go — -0.6 in glj according to WMAP7 [2(J. For smaller 
values of q, i.e. closer to —1, the inner region (red) would increase 
its area and the outer region (blue) would become smaller. For the 
value q = — 1, the blue region vanishes and there will be only red 
zone which corresponds, according to the condition q < —1 to all 
timclike directions having a positive contribution to the RE from 
the space-time geometry. 



Let us recall that a value of the deceleration parameter 
q equals to —1 is obtained for a(t) oc e Hot which is the 
solution in Robertson- Walker cosmological model of the 
Einstein's equations with cosmological constant A = 3Hq 
in vacuum. This solution is usually referred as the de 
Sitter expanding universe. This means that a value of 
q lower than —1 implies an expansion rate higher than 
the rate of expansion of a de Sitter solution. Therefore, 
it seems very unlikely to have q < — 1 in any epoch of 
the evolution of the Universe. Let us recall that even in 
the inflation phase the rate of expansion is expected to 
be exponential and therefore conflation = —1- Thus, let 
us conclude this section by remarking that for timelike 
vectors sufficiently close to the light cone one obtains 
< in any epoch of the evolution of the Universe. 

IV. CONCLUSIONS 

In this investigation we have studied the RE and its 
implications on the evolution of geodesies. We addressed 



the kinematics of a congruence of timelike geodesies with- 
out vorticity in a Robertson- Walker cosmological model 
with no underlying gravitational theory assumed a pri- 
ori. From a general point of view, we have established 
the difference between the implications of a cosmic accel- 
eration a > and a positive contribution to the RE from 
the space-time geometry, i.e. = —R a b£, a £, b > 0. 

We have then considered the fundamental congruence 
generated by the vector field £rw- We proved that the 
tendency to focus or diverge for the fundamental con- 
gruence depends of the fact that q is either bigger or 
smaller than — 1 respectively. Given the fact that an 
approximate value for the deceleration parameter today 
<7o ~ —0.6 can be found either by assuming GR or just 
by using Supernovae la data in theory-independent cal- 
culations, one concludes that although the fundamental 
congruence is expanding because of 9rw — 3d /a > 0, it 
has a tendency to focus since dO^w/dt < 0. 

Moreover, we have examined the congruence of time- 
like geodesies generated by a vector field orthogonal to a 
one-parameter e dependent spacelike hypersurface. The 
spatial coordinate velocity /3 = e/a 2 (t) represented - 
at each point of the hypersurface - the inclination of 
the vector field with respect to the preferred direction 
of Robertson- Walker models £^ w . In the case of the 
slice of a vector field to be normalized and orthogonal 
to the hypersurface, it is possible to determine the con- 
dition for the deceleration parameter q which ensures 
the contribution of the space-time geometry to the RE, 
Ai^a = —R a b£, a £, b to be positive. The aforementioned 
condition depends both upon the spatial coordinate ve- 
locity (3 and the scale factor a(t) but only through the 
physical spatial velocity a/3. In congruences without vor- 
ticity, if > docs not hold all the contributions to 
the RE are non-positive since w a b = and the remain- 
ing terms are non-positive as we discussed. We found 
that this condition is not satisfied in the present Uni- 
verse provided that a/3 > 0.8. This means that con- 
gruences of geodesies without vorticity with its tangent 
vector fulfilling af3 > 0.8 will not tend to diverge. In the 
Robertson- Walker preferred framework this situation oc- 
curs for geodesies moving near the light cones. These 
congruences will tend to focus provided that the other 
terms of the RE are negative. Besides, a negative contri- 
bution from space-time geometry to the RE would mean 
that the attractive character of gravity is restored at high 
speeds for — 1 < q < 0. The term high speeds must be 
understood with respect to the preferred framework of 
a Robertson- Walker cosmological model which is usually 
identified with the CMB because of its high isotropy. 

In the previous sense, systems may decouple from the 
cosmological accelerated expansion. This result can help 
to understand the structure formation of system moving 
fast, that is favored with respect to slow systems. Even 
taking into account the present value of the deceleration 
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parameter qo, a/3 > 0.8 gives a negative contribution to 
the RE from space-time geometry. This means that in 
order to get a negative contribution to RE from the space- 
time geometry, it is necessary to move with a velocity 
bigger than 0.8 c with respect to the preferred framework 
of Robertson- Walker space-times. Although this value is 
quite large, the deceleration parameter q was smaller in 
the past and consequently the velocity respect the CMB 
necessary to accomplish a negative contribution to the 
RE from space-time geometry was also smaller. 

To conclude, we would like to highlight the fact that 
the results obtained in this investigation provide us with 
a mechanism for a gravitational system to decouple from 
the accelerated expansion of the Universe. This mech- 
anism does not need to appeal to the gravitational in- 
teraction of the system itself. Although the decoupling 
of astrophysical systems from the cosmic acceleration is 
usually carried out by the stronger effect of their own 
gravitational forces, our investigation may help to un- 
derstand how geodesic motion in an expanding Universe 
influences this decoupling. 
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